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Abstract 

Requirements for correlation measurements in high-multiplicity events are discussed. 
Attention is focussed on detection of so-called hot spots, two-particle rapidity 
correlations, two-particle momentum correlations (for quantum interferometry) and 
higher-order correlations. The signal-to-noise ratio may become large in the 
high-multiplicity limit, allowing meaningful single-event measurements, only if the 
correlations are due to collective behavior. 
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1. Introduction 

In the next ten years, ultra-relativistic heavy ion collisions are planned at 
Brookhaven's Relativistic Heavy Ion Collider (RHIC) and CERN's Large Hadronic 
Collider (LHC). In both cases, experimenters expect to see multiplicities in excess of 
1000 particles per unit rapidity. As a result, there has been a lot of speculation about 
single-event fluctuation measurements [|T]-[J. In this paper, I assess requirements for 
various single-event measurements, and calculate the number of events needed for 
useful measurements in cases where single-event analyses are not meaningful. 

I begin with a discussion of searches for so-called hot spots (regions of unusually 
high particle density) in Section ||]. In Section [| I discuss measurements of two-particle 
rapidity correlation functions, assuming that high-multiplicity events are independent 
superpositions of lower-multiplicity events. In Section £|, I discuss measurements of 
two-particle momentum correlation functions that are commonly constructed to 
measure collision volumes using quantum interferometry. In Section [5], I compare 
results from two-particle correlation functions with those from higher-order correlation 
functions. Finally, I summarize the results in Section || 



2. Looking for hot spots 

One common suggestion is that it may be possible to search for hot spots, or regions 
with unusually large numbers of pions. The basic motivation for these searches is 
simple: any process that creates a lot of entropy in a small rapidity bin is of interest. 
Thus, hot spots are commonly thought of as possible signals for interesting phenomena. 

For defmiteness, suppose that searches are made within a rapidity window of size 
AY, using events with N particles in this window. If particles are randomly distributed 
with a flat rapidity distribution, the mean number of particles in a bin of size Sy is 

where dN/dy = N/AY. The standard deviation is 

a n = - n 2 = n{l-5y/AY) « n, (2) 

if the window is large (Ay ^> 5y). 

The central limit theorem applies in the limit iV — > oo, so particle number 
fluctuations have a gaussian distribution. The probability that a given bin contains 
more than n + Sn particles is 

if Sn 2 ^> n. To achieve success in a hot spot search, the hot spots must be present 
significantly more often than in a random distribution: 

f(5n) > V(5n/n^ 2 ), (4) 
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where f(Sn) is the probability that a given bin contains a hot spot with at least 5n 
excess particles. 

For definiteness, I assume that a useful result must find hot spots at least ten times 
as often as expected. In this case, hot spots that occur in 10% of the bins must produce 
at least 2An 1 / 2 excess particles. For RHIC and LHC events, particles are produced 
over about ten units of rapidity, so hot spots that occur in more than 10% of bins will 
be seen more often than once per event, and are thus not very useful as triggers for 
interesting events. The larger dN/dy is, the larger the hot spots must be before they 
can be separated from the background fluctuations, so it is likely that hot spot searches 
will be most profitable in events of high energy but relatively low multiplicity. 

For RHIC and LHC experiments, rapidity densities in excess of dN/dy = 1000 are 
expected. A hot spot that produces high energy pions (p > m in the hot spot rest 
frame) isotropically has a width of at least one unit of rapidity, so I take 5y — 1. In 
this case, if hot spots occur in 10% of the bins (about one hot spot per event), they 
must produce 75 excess charged pions to be clearly useful as a trigger. If they occur in 
1% of the bins (about one per 10 events), they must produce 100 excess pions. 

One possible mechanism for visible hot spots is the production of bubbles of 
disordered chiral condensate These bubbles were proposed as an explanation for 
the so-called Centauro events observed in cosmic ray studies, in which many charged 
pions were produced with very few neutral particles, in contrast to typical nuclear and 
high-energy processes that produce equal numbers of 7r + , 7r~, and 7r° mesons. A 
chirally-disordered bubble that produces iV ± charged pions, with no 7r° mesons, yields 
N ± /3 excess charged pions within one unit of rapidity of the bubble. If these bubbles 
are produced in less than 10% of the bins, they will be clearly visible if N^ 1 > 225, but 
this is unlikely as the expected value is N ± « 20 Q. Thus, hot spot searches at RHIC 
and LHC probably will not yield evidence for bubbles of disordered chiral condensate. 

Quark-gluon plasma (QGP) droplets are typically too small to be hot spot 
candidates in high-multiplicity events. A typical mean radius for a QGP droplet is 1 
fm, at which size the droplet should produce about 18 charged pions. There is no 
proposed mechanism that would produce hot spots large enough to be clearly 
observable at RHIC or LHC; however, if they are seen, the lack of theoretical 
prediction will make them even more interesting than if they were expected. 



3. Two-particle rapidity correlations 

In this section, I discuss measurements of two rapidity correlation functions: the 
standard two-particle correlation function, R2, and the simplest split-bin correlation 
function, S^. The standard two-particle rapidity correlation function is |5|. |(J 

where AY is some large rapidity separation that is used as a reference. I assume at first 
that p is flat, and discuss the effect of corrections for non-flat distributions afterwards. 
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Suppose that I take some arbitrary model of particle production and analyze a single 
event. Consider events with N particles in a rapidity window of width AY, where the 
rapidity distribution for any given particle is p(y) = 1/AY. Let there be N c correlated 
pairs, where typically iV c <C N(N — 1), and the rapidity distribution for correlated 
pairs is q(yi — t/2)/AF, where q is some arbitrary function. For simplicity, I consider 
only events with exactly iV particles in the rapidity window; generalization to events 
with differing multiplicities is straightforward 0. 

I do not assume that all correlations are pair-wise, but I neglect higher-order 
correlations for the moment. For example, it is possible that the pair-wise correlations 
result from interactions of large numbers of particles. Even in this case, however, 
correlation functions are dominated by pair-wise correlations unless the interactions 
involve almost all of the particles. I discuss this in more detail in Section [|. 

Finally, I assume that a superposition of a independent events with n particles each 
is equivalent to a single event with an particles. This is equivalent to assuming 
independent nucleon collisions, or independent parton collisions. In this case I must 
have N c = kN, where k is some unknown proportionality constant. [For example, if I 
combine two events I double both iV and N c , as pairs of particles from different events 
are clearly uncorrelated.] 

I can immediately write down the two-particle density, 

^y^ WN - 1) - N ^ Y ^ Yq(vi - y > ) . (6) 

Using Eq. (||), I obtain the two-particle correlation function, 

R 2 (y; AY) = 1 + -Jf— [q(y) - q(AY)} . (7) 

Here (and for the remainder of this paper) I drop corrections of order 1/N and 1/AY 
unless otherwise specified, as I am primarily interested in the analysis of 
high-multiplicity, high-energy events. If the mean separation for a correlated pair is y*, 
then typically q(0) ~ while q(AY) — > for large Ay, so the maximum value of 
R 2 is roughly 

k 

I obtain a lower limit for the error in a measurement of R 2 by calculating the 
expected fluctuations in the absence of correlations. Consider an experimental 
measurement: 

Nev 

^2 7ii(0) n t (y) 

MvAY) = j^ • (9) 

£ m(0) m(AY) 
i=i 

Here rii(z) is the number of particles in the i-th event with rapidities between z — 5y/2 
and z + 5y/2 (Sy is thus the experimental bin size), and N ev is the number of events 
used in the measurement. 
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Assuming that Sy < y, so that the bins do not overlap, and that there are no 
correlations, 

N(N-l) 5y 2 



(n t (0)n t (y)} 
independent of y, so clearly (R 2 ) 



AY 2 

1. The standard deviation is 



(10) 



N c 



E ^ 2 (o)^) 



i=l 



+ 



£n 2 (oK(Ar) 

i=i 



E MO) ni(y) 

Li=l 



.t=i 



(11) 



In the absence of correlations, the standard deviation is 

(-8JV + 12) Sy 2 + 4(JV - 2)AY5y + 2AY 2 



N ev N(N - l)Sy 2 



(12) 



For a good measurement, the experimental bin size must be much smaller than the 
total rapidity window used, so Sy -C AY. As the purpose of this paper is to consider 
measurements in high-multiplicity events, I take the limit iV —>■ oo, and obtain 



4 AY 



1 + O 



AY 



N ev N5y 

Thus, the error in the measurement is 

2/y/N ev (dN/dy) Sy 



(13) 



(14) 



where dN/ dy = N/ AY is the rapidity density. 

It is possible, if the fluctuations in the system are large, that the actual error is 
larger than given by (|T4"D. If the measured standard deviation is smaller, however, then 
the value of eq. (|I~4"D is probably better to use, as this represents the error in measuring 
uncorrelated events. If the measured fluctuations are anomalously small, then the 
system is probably strongly correlated, so if i? 2 — 1 is not significantly different from 
zero it is probably best to look for another correlation function that reflects these 
strong correlations. 

Combining eqs. (§) and (|i~4|), I obtain the signal-to-noise ratio, 



[s/n) 



k 



R 



( N ev Sy \ 
2y* \dN/dy) 



1/2 



(15) 



If all particles are produced in clusters containing n c ^> 1 particles, and all particles in 
a given cluster are pair-wise correlated, then k « n c . [This is trivial - every particle is 
produced with n c associated (correlated) particles, so there are n c correlated pairs per 
particle.] Finally, for a passable measurement Sy < y*, so the best possible 
signal-to-noise ratio is 



{s/n)! 



N P 



nc ( 

2 {(dN/dy) y* 



1/2 



(16) 
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I use eq. fllq) to estimate how many events I need in order to measure R 2 well, as a 
function of the cluster size. For most clusters, y* < 1.3, so to measure R 2 (at its peak) 
with better than 4<r accuracy I need 



1/2 



2 V 1.3 dN/dy 



>4, 



(17) 



or 



N ev > 83 (dN/dy) / n, 



(18) 



Here I use y* = 1.3; this is obtained for the most energetic clusters, and gives the most 
pessimistic estimates of s/n. The cluster size seen in nuclear collisions || |9|] at 200 
GeV is approximately ten charged particles [TIJ. If this persists up to RHIC and LHC 
collision energies, where the charged particle multiplicity dN c h/dy ~ 1000, then 
approximately 830 events will be needed to obtain a good measurement of the peak 
value of R 2 . 

This is, of course, a naive theorist's estimate, leaving out any possible experimental 
difficulties, and applies only to a measurement of the amplitude of R 2 . If I want to 
measure the shape of R 2 reasonably well, I should really require that 5y < y*/10, in 
which case I find that I need 8,300 events. However, eq. (|1^) does illuminate the 
difficulty of measuring correlation functions in high-multiplicity events: if cluster (or 
source) sizes are independent of dN/dy, then s/n decreases with increasing multiplicity, 
and accurate measurement becomes increasingly difficult. Also, eq. (|1(J) justifies a 
posteriori the neglect of correlations when calculating the measurement error, as for 
one event the correlations (the signal) are much smaller than the statistical fluctuations 
(the noise). 

One could argue that the measurement I have outlined for R 2 does not efficiently use 
the available statistics. As a response, I construct split-bin correlation functions 
(SBCFs) [PH] , in order to use the available statistics with maximal efficiency. The 
simplest second-order SBCF is 



S 2 (6y,AY) 



AY l^ y fU-imSy fjSy 

AF £ / d Vl dy 2 pW( yi ,y 2 ) 

by \ dy x / dy 2 p (2) (yi,y 2 ) 
Jo Jay/2 



(19) 



Here AY/5y must be an integer; taking 5y = AY/2 1 for i — 0, 1, 2, . . . uses all of the 
two-particle phase space without re-using any pairs of particles. 
Under the assumptions that I used to calculate R 2 , 



S 2 (5y;AY) = 1 



k 



dN/dy 



[g(5y) - g(AY)) , 



(20) 



where 



9W 



z/2 rz 

dxxq(x) + / dx (z — x) q(x) 

ii Jz/2 



(21) 
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For Sy < y*, if q(z) is linear in z, 



g(5y) = q(5y/2), (22) 

while for quadratic q(z), 

g(5y) = q(5y/^24/7). (23) 

Eq. ( |23"1) is very close to eq. (|2~2"D , so these relations are insensitive to the shape of q and 
are thus robust. For AY S> 2y*, 

9 (AY) = ^ - 0, (24) 

so for Ay > 2y* and 5y < y*, 

S 2 (5y-AY) « i? 2 (^/2), (25) 
independent of Ay. The maximum value of 5*2 (Ay) is then 



S™ ax {AY) = 1 + 



*2 



AF 2 



A- 



(diV/dy) j/* 



(26) 



I calculate the error in the same manner as before, assuming that there are no 
correlations: 

4 

as ~ N ev N 
Thus, the error in the measurement is 




(27) 



e s = 2/^N ev (dN/dy) AY, (28) 
and the signal-to-noise ratio is 

/ / \ k (N ev AY\ 1 ' 2 , s 

for AY > 2y*. Setting k = 10, y* = 1.3, AF = 10 and dN/dy = 1000 as previously, I 
need approximately 110 events to measure S 2 to 4cr accuracy, as compared with 830 
events to determine the peak value of R 2 - 

The reader should note that this is an estimate of the minimum requirements for a 
measurement of S 2 - Applying eq. (p9]) to a sample of 92 central O+Em events at 200 



GeV, with y* = 1.3, dN/dy w 40, and AY = 4, I obtain (s/n)s ~ 12. If I use eq. fl2§) 
to estimate the corrections for the finite value of Ay, I obtain (s/n)s ~ 7. This is in 
reasonable agreement with the observed value (s/n)s ~ 2 — 3 [Q. Most of the 
difference comes from the crude approximation used for ffi^, as an exact calculation 
T3fl shows that eq. (@) overestimates the signal by a factor of two. Thus, the estimate 



of s/n is within a factor of 3, while the number of events needed for a good 
measurement is a factor of about 10 more than estimated. 
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The measurement of S2 discussed above would give the shape of R2 with points that 
are approximately 0.7 apart on a logarithmic scale, as I keep changing the bin size by 
factors of two. The previously discussed measurement of R 2 gave points approximately 
0.1 apart (spacing was y*/10). Duplication of this measurement using S2 would involve 
seven independent measurements of S2, requiring approximately 800 events, as opposed 
to the 8,300 required for the R 2 measurement. Thus, if the same number of events are 
used in each case, a measurement of S2 has a bit less than one-third of the statistical 
noise of the corresponding measurement of i?2- 

If the single-particle distribution is not flat, the correlation functions, R 2 and S 2 , 
should be modified. The most useful two-particle correlation function is 

pW(0,y)p(AY) 
R ^ AY) = pV(0,AY)p(yy (30) 

„ p (2) (o,y)(iv) 2 (3l) 



p(0)p(y) (N(N-l)) 

These rapidity correlation functions are unity whenever particles are distributed 
randomly in rapidity according to the single-particle distribution, independent of the 
multiplicity distribution. I outline the revised calculations for R 2 below. 

Suppose that the distribution of cluster centers is n ((yi + y 2 )/2), instead of simply 
1/ AY as assumed so far. If all particles come from clusters, then the shape of the 
single-particle distribution is not flat. The probability that a particle is found between 
y and y + dy is then p(y)dy, where 

p(y) = = / dyMvc) q(2y c - 2y) » *(y), (32) 

if n(y) varies slowly compared to q [q{2y c — 2y) occurs because the second particle of a 
pair has rapidity y = 2y c — y\. I then obtain 

p (2 \y u y2) = [N(N - 1) - N c }p( yi )p(y2) + N c n (^^) ^ " ^ 2 )' ( 33 ) 
and consequently 

R '^ AY) ^ + W^' (34) 

r i 2 

assuming that dN/dy varies slowly so [dN/dy\ ] [dN/dy\ y ] ~ dN/dy\ y /2 ■ 

The only difference from R 2 is that dN/dy is now taken at rapidity y/2. The error 
will change slightly, as the error in determining the correction factor should now be 
added to the previous error. The error in determining p(y) is 



e p /p = l/^/N ev (dN/dy) y 5y, (35) 

so the error in R' 2 is naively 



e R = y/6/N ev {dN/dy)Sy, (36) 
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assuming that the distribution is approximately flat and that the errors are 
uncorrelated. This will lower the previous values of (s/n)#, and raise the numbers of 
events required for a good measurement by about 50%. 

The most useful version of 5*2 has only a simple additive correction [12]: 

S' 2 (5y; AY) = S 2 (5y; AY) - Z 2 (5y; AY), (37) 

where 



AY l & y r(j-l/2)Sy rjSy 

AY J2 dyip(yi) / dy 2 p(y 2 ) 

MSy;AY) = 3 1. (38) 

5y / dyip(y!) / dy 2 p{y 2 ) 
Jo Jay/2 

For the case where the particles are centered in the window, 
Jo^ Y ^ 2 dyp(y) = Jay/2 dyp(y), the correction simplifies to 

4 Ay / 5 f r(j-l/2)Sy rjSy 

Z 2 (6y;AY) = — — £ / dy l Sf(y 1 ) dy 2 6f(y 2 ). (39) 

AY5y fr[ J(j-l)Sy J(y-l/2)Sy 

Here 5/ = pAF — 1 is the fractional difference from the mean value. As by — > 0, 
^2 ~~ ^ <^/ 2 from below, so if the furthest bins are within fraction / of the mean density 
then S 2 < f 2 . The error in the correction factor, S 2 , is approximately the same size as 
the error in S 2 , so the experimental errors will be somewhat larger than for the 
uncorrected case, just as for R 2 . 



4. Quantum interferometry 

Experimenters are also very interested in using Hanbury/Brown-Twiss interferometry 
to determine the geometry of the collision region in high-multiplicity events Jl4|, [15|]. 
This technique uses species-dependent momentum correlation functions, for which 
particle identification is required. I briefly discuss statistical aspects of the two most 
common methods for identical particle interferometry, neglecting any difficulties due to 
errors in the identification of particles or other technical problems. 

Both techniques usually use pions that are identified in a spectrometer covering A$ 
radians of azimuthal angle and AY units of pseudo-rapidity. For pions, 
pseudo-rapidity and rapidity are almost identical, so I do not differentiate between 
them in this paper. The correlation function used is 



J A$ 6^02 J Ay dy x dy 2 J dp^dpT^ 5[q 2 + (gi - p 2 f 
J A * d(j)idcj) 2 Jo AY dy x dy 2 J dp T ,idp T , 2 puc 8[q 2 + (pi - p 2 f 



(2) 

where 5 is the Dirac 5-function. Here p id is the distribution for two identical particles, 
and p^ 2 } is an uncorrelated two-particle distribution; the two methods differ only in the 
prescriptions used to construct p^. I use standard high energy units with h — c — 1. 
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In the first method, referred to as event mixing, an uncorrected distribution is 
produced by combining particles (of the same species and charge as those used in p id ) 
from different events. This is usually done by constructing simulated events using 
particles from measured events, while ensuring that no two particles come from the 
same events. This process is very computationally intensive for high-multiplicity 
events, as it is necessary to keep track of the event from which each particle is taken 
and check all selections to ensure that they don't come from an event that was used 
earlier. 

Considerable computational difficulty is removed by constructing pP} by convoluting 
the single-particle distribution obtained by averaging over all events. Using a sample of 
N ev events, each with N particles in the spectrometer, this second procedure yields 

AT _ 1 „(2) _ (2) 

P (2) (Pl,P2) = p(pi) P (P2) = A»S + ^ N J UC , (41) 

for all pi 7^ p2, where p^} is the value obtained by taking all convolutions with no two 
particles from the same event. For N ev 3> 1, the difference between the proposed 
procedure and the usual one is small. This change in procedure is even more important 
for measuring higher-order correlations, as the construction of p uc by the usual event 
mixing quickly becomes computationally prohibitive, while the simple convolution of p 
is almost always feasible. 

In the second method, referred to as charge mixing, is constructed using 
particles that are identical except for charge. For example, 

A>2(pi,ft) = P + (Pi)P~(P2), (42) 

where p ± is the n^ density, is often used for pion interferometry. As charge mixing 
does not require information from more than one event, it is the preferred 
normalization technique for analysis of single events. 

Consider the analysis of a single event. The mean number of ^ mesons seen in a 
spectrometer covering AY units of rapidity and A$ radians of azimuthal angle is 

N± = (dN/dy) AY A<P ^ 
An 

if particles are spread randomly in y and <fi with a uniform distribution. I imagine an 
ideal spectrometer that detects all pions passing through it, independent of transverse 
momentum, pt- For simplicity, I assume a thermal-type px distribution: 

V( PT ) = ^ e ~ 2PT/{PT) > ( 44 ) 

where V(pt) dpr is the probability that a given pion has transverse momentum 
between px and px + dpr, and (pr) is the mean pt- 

For a spectrometer large enough that edge effects are unimportant 
(AY, A<3> 3> q/ipr)), the standard deviation of C 2 (g) is 

167T f (I 2 } \ 

° C ~ (dN/dy) AY A$ \ (J) 2 + 2 (dN/dy) qSq(I)j' { } 
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in the high-multiplicity limit, where 



/•CO roc roc rlix 

(I n )= dp T7l V(pT,i) / dpT >2 V{jpr,2) dy 

JO JO J-oo Jo 



5 (q 2 + 2m 2 + 2ptaPt,2 cos (f> — 2 ^ (px t i + m2 )(p'r ) 2 + m2 ) cosh, y) . (46) 

Here (/ 9 9 =0 dq 2 ) n (I n ) / (AY A$) n is the probability that n given particles each have 
momenta within q' of a particle at y = <ft = 0, with transverse momentum distribution 
V(pt,i)- The primary interest is in small q, and m 2 -C (pt) 2 , so I evaluate (I n ) for 
g 2 = m 2 = 0, obtaining 

and thus 

!67r [4 y^Fm(p T ) 1 

C '"' (dN/dy) AY A$ \ 7r [r(l/4)] 2 (diV/cfy) g5g J ' 1 1 

The first term of eq. ( f49|) dominates as long as 

q 2 » ^ (PT) , (50) 

[r(l/4)] 2 (dN/dy)' 

where I have taken Sq = q for the smallest bin. Using dN/dy = 1000 and (p?) = 500 
MeV, I find that the first term dominates as long as q 2 30 MeV 2 , which is true for 
all practical measurements. Thus, the expected measurement error is 

e c « , 4 (51) 
'(dN/dy) AY A$ 

independent of the momentum or the bin size! As the signal is approximately unity for 
q = 0, a 4<7 determination is possible with a single event for any spectrometer with 
AFA$ > 256 (dN/dy)- 1 . 

Quantum interference also produces two-particle rapidity correlations. These 
correlations are due to collective effects, so they have a different multiplicity 
dependence than the two-particle correlations discussed in the previous section. The 
momentum scale for quantum interference q* ~ 1/r, where r is the size of the system, 
so the rapidity scale is 

a* 4 fm 

V* « 7^ « — , (52) 
(Pt) r 

while the number of correlated pairs per particle is 

k ^ g* 2 (I) (dN/dy) ^ 0.5(tZAVtfy)fm 2 
47T r 2 
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For simplicity, I assume that pairs are correlated if the momentum difference is less 
than q*, and uncorrelated otherwise. Using eqs. ( |20| ) and fl24|), I obtain 



for Ay ^> Sy 3> lfm/r, while eq. (|26|) gives the maximum value, 

be , max ^ i + L25ftn_ 
r 

The two-particle correlation due to clusters is visible above S^ 6 as long as 
dN/dy < 1.25n c y*r 3 /fm 3 . Single-event measurements of S^ 6 are in principle possible 
if (dN/dy) AY > 40r 2 /fm 2 . Recent data [TO indicates that r = 3 — 4 fm for central Si 
collisions at 14.6 GeV/nucleon (approximately the Si radius); using r = 7 fm for U+U 
collisions at RHIC and LHC, a single-event measurement of 5*2 may be possible. 
However, such a measurement is very difficult technically, requiring rapidity resolution 
to y* » 0.06. 



Higher— order correlation functions 



It is also possible that higher-order correlation functions might give better results 
(or more interesting results) than two-particle correlation functions. Higher-order 
correlation functions can in principle be used to determine three-body and 
higher-order interactions, and many experimenters have tried to use them for this 



purpose, although without much success [[17], [18| . Alternatively, measuring 
higher-order correlation functions might provide a more accurate determination of the 
two-particle correlation function than can be obtained from a direct measurement. 

To test these hypotheses, I construct scaled factorial moments (SFMs) ||19|| , that also 
use the data more efficiently (although they re-use pairs of particles). For pedagogical 
purposes, I consider only the so-called exclusive SFMs, 



Fi(5y;AY) 



(AYY' 1A ^ y Ju-m 

\5y) h N ■ 



dyi ■•■ dyip i% \y u ...,yi) 



i=i 



(N-i + 1) 



(56) 



Just as for S2, AY/5y must be an integer; however, any two values of 5y use some 
common phase space, so it is impossible to construct SFMs without over-using the 
available phase space. 

To examine the feasibility of studying higher-order correlations, I extend my 
previous approximations to include three-body correlations, assuming N t correlated 
triplets with the boost-invariant distribution q t (yi — 2/2, Hi — 2/3)/ 'AY '. For AY ^> y*, 
the maximum value of F,- is 



1 + 



1) 



k 



(dN/dy)y* 



+ 



ijj, - l)(i 



k t 



J [{dN/dy)y* 



2 ■ 



(57) 
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where k t = N t /N. For cluster decay, k t ps n\ for n c 1, so 



(»(»-!)(» -2) 
6 



Av 



2fz - 2) 



J [(dN/dy)yf - 1) 



(58) 



is the three-particle contribution to Fi, while Fi — 1 is the two-particle contribution. 

The observed three-particle correlation decreases with increasing multiplicity as 
(dN/dy)~ 2 , even faster than the two-particle correlation. Equation (|58|) is apparently 
very general, as it holds for any values of n c and y*. Given a distribution of values, 
corrections of order unity are likely; however, it is still probable that the difficulty of 
extracting the three-particle correlation increases very fast as the two-particle 
correlation function decreases. This is, in my opinion, the most likely reason for the 



failure of experimenters to extract significant three-particle correlations fli7| , |18f from 
their data, as virtually all two-particle correlations are approximately 10% or smaller. 

The second hypothesis, that higher-order correlation functions might provide a more 
accurate determination of the two-particle correlations than direct measurement, seems 
to be true. Under the assumptions that I use to calculate R 2 , 



Fi(6y) = 1 + 



in 



k 



dN/dy 



[h(5y) - h(AY)} 



where 



2 r z 

h(z) — — I dx(z — x)q[ 
z l Jo 



x) « q{z/3). 



I calculate the error in the same manner as before, assuming that there are no 
correlations: 



1 + O 



Sy 
AY 



(59) 



(60) 



(61) 



The statistical error in Fi is proportional to i, while the signal is proportional to 
i(i — 1), so s/n is proportional to i — 1 and thus improves with increasing i. However, 
the above arguments are valid only when (dN/dy) 5y 3> i, in which case most bins 
contain i particles, and even then apply only to statistical noise. In most cases, there is 
also noise from undesired correlations produced by the detector, and this systematic 
noise is proportional to i(i — l)/2. It is thus possible that s/n is approximately 
independent of i. 

It appears that measuring higher-order correlation functions is the most accurate 
way to determine two-particle correlations. Because SFMs re-use data, the cleanest 
approach is probably to use higher-order SBCFs [|llj] : 



Si(6y;AY) 



AY/Sy i y_ in _ 1)/i]Sy 

ay £ n / W'W-.w) 

j=l n=l J [i-n/i\8y 



pnAY/i 

Sy 1 ' 1 II / dy n p^\ yi 

„ = 1 J{n-l)AY/i 



(62) 



Vi 



Maximal use of the data without re-use is simple for S 2 , but not easily achieved for 
high-order SBCFs. As a result, it may be preferable to use some different 
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generalization of S2 for maximal efficiency. Higher-order correlation functions might 
also give better results for quantum interference measurements than the commonly 
used two-particle correlation functions. 

6. Conclusions 

I have discussed four types of correlation measurements: hot spot searches, 
two-particle rapidity correlations, two-particle momentum correlations (for quantum 
interferometry) , and higher-order correlation functions. Hot spot searches are most 
likely to be profitable in events of high energy but relatively low multiplicity. 
Two-particle rapidity correlations are most easily measured in events of relatively low 
multiplicity, if high-multiplicity events are just superpositions of lower-multiplicity 
events. A good measurement of S2 at RHIC or LHC will require at least 800 events. 

Single-event measurement of two-particle momentum correlations due to quantum 
interference is possible in high-multiplicity events with spectrometer coverage 
AFA$ > 256 (dN/dy)-\ which should be easily attainable at RHIC and LHC 
Rapidity correlations due to quantum interference are in principle measurable in single 
events at RHIC and LHC, but such measurements would be very difficult technically. 
Measuring higher-order correlation functions in high-multiplicity events gives little 
information about three-body and higher-order correlations. However, measuring 
higher-order correlation functions can give a more accurate determination of 
two-particle correlations than direct measurement of two-particle correlation functions. 
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